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Abstract 

We investigate the transport of Brownian particles in a two-dimensional potential under the 
action of a uniform external force. The potential is periodic in one direction and confines the 
particle to a narrow channel of varying cross-section in the other direction. We apply the standard 
long-wave asymptotic analysis in the narrow dimension and show that the leading order term is 
equivalent to that obtained previously from a direct extension of the Fick-Jacobs approximation. 
We also show that the confining potential has similar effects on the transport of Brownian particles 
to those induced by a solid channel. Finally, we compare the analytical results with Brownian 
dynamics simulations in the case of a sinusoidal variation of the width of the parabolic potential 
in the cross-section. We obtain excellent agreement for the marginal probability distribution, the 
average velocity of the Brownian particles and the asymptotic dispersion coefficient, over a wide 
range of Peclet numbers. 
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INTRODUCTION 



Recent progress in microfluidic devices has led to the development of novel separation 
strategies that take advantage of the unprecedented control on the geometry and chemistry 
of the stationary phase at scales that are comparable to the size of the transported species 
. A fundamental problem that is at the core of several of the proposed separation 
devices is the transport of Brownian particles through entropy barriers. A representative 
example is the motion of a suspended particle in a channel with periodically varying cross 
section. The purely-diffusive transport in the absence of an external force has been studied 
extensively [4, y, [7], and a well-known approach is to reduce the dimensionality of the 
problem via the Fick- Jacobs approximation. In this approximation the motion in the cross 
section is reduced to an entropic contribution to the longitudinal transport 8]. The case of 
biased diffusion in the presence of a driving force has also received considerable attention due 
to its relevance to separation devices. Hanggi and coworkers examined the validity of a direct 
extension of the Fick- Jacobs equation to describe the biased motion of a Brownian particle 
in a narrow channel of varying cross section. The authors showed that using an ad-hoc 
position-dependent diffusivity provides a good approximation, for Peclet numbers below 



also investigated the biased motion 



a critical value 0, [h]]. Dorfman and coworkers 
of a Brownian particle in a periodic channel. Specifically, they performed an asymptotic 
perturbation analysis to obtain simple expressions for the macroscopic transport coefficients 
that remain valid at relatively large Peclet numbers. 

Here, we also consider the transport of Brownian particles confined to a channel of period- 
ically varying cross section but, in the present case, the confinement is induced by a potential 
energy landscape and not the solid boundaries of a channel. This type of spatial confinement 
occurs, for example, in the transport of suspended particles in microfluidic channels. De- 
pending on the density of the particles, the van der Waals and electrostatic forces between 
the particle and the channel walls, and the chemical composition of the media, the particles 
could become confined to the secondary minima of the particle-wall interaction potential in 



the vertical direction 



12] . In that case, the presence of a periodic pattern on the bottom 



wall, such as that created by the deposition of thin metal stripes perpendicular to the flow, 
would lead to the confinement of the Brownian particles to a periodically varying channel 
parallel to the wall, analogous to the one considered here. We shall show that the effect on 
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the average velocity of the particles (and the probability distribution in general) induced by 
this type of soft confinement of the suspended particles to an energy minima is the same as 
that caused by geometric confinement between solid walls. Understanding the effect that 
this type of confinement has on the transport of suspended particles is particularly impor- 
tant for the development of recently proposed separation techniques in microfluidic devices 
that are based on partitioning jl3| . 



TRANSPORT OF BROWNIAN PARTICLES IN A CONFINING PERIODIC PO- 
TENTIAL 



Let us consider the transport of Brownian particles in a potential that is periodic in 
the x-direction, V(x = x ,z) = V(x = x + L,z), and that confines the particles in the 
z-direction, that is V(x, z) — ► +00 for z —>■ ±00. In the limit of negligible inertia effects 
the motion of the particles is described by the Smoluchowski equation for the probability 
density P(x, z, t), 

^ + V-3 = 5(x,z)5(t). (1) 
The probability flux, J(x,z,t), is given by 

J = I ( FP - - kB T 9 ?) U U-^-P- k B T d -f] k, (2) 



77 \ dx dx J r) \ dz dz 

where F is a uniform external force in the x-direction, 7] is the viscous friction coefficient, 
and we have used the Stokes-Einstein equation to write the diffusion coefficient as a function 
of 77, D = ksT/rj. 

In order to obtain the asymptotic distribution of particles within a single period of the 
potential we first introduce the reduced probability density and the reduced probability 
current (see Refs. Q, Q or the analogous approach presented in Ref. 16]), 

+00 

P(x,z,t)= P(x + n x L,z,t), (3) 

n x =—oo 
+00 

3(x,t) — J(x + n x L,z,t). (4) 

n x =—oo 

The reduced probability is obtained by solving the Smoluchowski equation with periodic 
boundary conditions in x. In particular, the long-time asymptotic probability density, 
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Poo{x, z) = lim^oo P(x, z, t), is governed by the equation (dropping the tilde), 



V- J™ = 0. 



(5) 



The far-field condition in z is a vanishingly small probability density and flux due to the 
confining potential, 



J, 



z 

oo 



l_(_dV 

1 o 1 oo 

OZ 



k B T 



dP^ 



-> 0. 



(6) 



7] \ OZ dz J z^ioo 

Finally, the reduced probability is obtained by imposing the periodic boundary conditions 
in x, 

P oo (x = 0,z) = P oo {x = L,z), (7) 



and the normalization condition, 

(Poo) 



L poo 

dx / PrscdZ 







1. 



(8) 



In the case of narrow or slender channels, we assume that the characteristic length scale 
in the cross-section, i. e. perpendicular to the channel centerline, is given by eL, where 
e 1 is the slenderness ratio. We can then write the governing equation and boundary 
conditions using the following dimensionless variables, x = x/L, z = z/(eL), V = V/(fcflT), 
where ks is the Boltzmann constant and T is the absolute temperature, and the re-scaled 
probability density P^ = eL^P^. For the sake of simplicity all bars are dropped after 
nondimensionalization. The governing equation for the reduced probability (Eq. then 
becomes: 



,d_ 

dx 



Pe 



dV_ 

dx 



Po 



dPoo 
dx 



d_ 

dz 



r oo 

dz 



dP 

ul oc 

dz 



0. 



(9) 



where the Peclet number is defined as Pe = FL/ksT. Simple inspection of this equation 
shows that the leading order approximation for e 1 is the local equilibrium in z, which 
corresponds to Jq(x, z) = 0, due to the no-flux far-field condition. In general, an accurate 
description of the long-time transport of Brownian particles can be obtained from the first 
two moments of the asymptotic probability, which describe the average velocity and the 
broadening of the distribution. Applying macrotransport theory we can write the average 
velocity in terms of the asymptotic probability distribution 161 ]. 



(v) 



dx 



Jtdz 



dx 



Pe 



dx J 



dPoo 
dx 



dz. 



(10) 
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The dispersion coefficient D*, can also be calculated from the asymptotic distribution via 
the so-called B— field, which is the solution of the following differential equation [h] ]. 



- (Poo — ) 

dz \ 00 dz J 



dV p _dP™ \ dB_ 

dz 00 dz I dz 



- (pJ—\ 

dx \ 00 dx J 



dx J 00 dx J dx 



e 2 Poo(v). 



The boundary conditions for the B— field are, 



dB 



0. 



dz z^±co 

B(x = l,z)-B(x = 0,z) = -1. 
Finally, the dispersion coefficient is given in terms of B(x, z) by 



D* 



dx 



dB 

dx 



1 (HE 

e 2 \ dz 



dz. 



(11) 



(12) 



(13) 



ASYMPTOTIC ANALYSIS IN THE NARROW CHANNEL APPROXIMATION 



We apply the standard long-wave asymptotic analysis to obtain an approximate solution 
to the problem described above. First, we propose a solution to the stationary probability 
distribution in the form of a regular perturbation expansion in the slenderness parameter, 



Poo(x, z) ~ p + e 2 pi + e i p 2 H . 

Analogously, we write a regular perturbation expansion for the probability flux, 

Joo(x, z) ~ J + e 2 Ji + e 4 J 2 H . 



(14) 



(15) 



Substituting these expansions into Eq. ([9]) it is straightforward to determine the governing 
equation for the leading order terms, 



d 



dV 



dp 



dz V dz ^° dz 



dz 



0. 



(16) 



The corresponding leading order boundary and normalization conditions, derived from 
Eqs. ([MED, are: 



Jq (x, z = ±00) 







(17) 



p (x = Q,z) = p (x = l,z) 
(P0> = 1. 



Then, integrating Eq. f|T6|) and taking into account the no-flux condition we obtain: 



p (x,z) = f (x)e- v ^ z \ 



(18) 



where fo(x) is an unknown function that is to be determined from the second order 0(e 2 ) 
balance of Eq. <M} , 



d ( dV dp\\ d 
dz \ dz ^ l dz J dx 



dx J 



Pe - — I po 



dpo 
dx 



(19) 



The corresponding boundary and normalization conditions are: 

J{ (x, z = ±00) = 0, 

pi (x = 0, z) = pi (x = l,z) 
(Pi) = 0. 

Substituting the solution obtained for p into Eq. (|T9"1) we obtain 
d ( dV dpi\ d 



(20) 



dz V dz ^ dz ) dx 



Pe/o - -1- 1 < 
ax 



V(x,z) 



(21) 



In steady state, the total flux in the x-direction is constant along the channel. Therefore, 
by integrating the equation above over the cross-section, and taking into account the zero 
flux condition in the z-direction, we obtain 



d_ 

dx 



Pe/o - £ ) /(*) 



0. 



where 



Kx) 



~ v ^dz. 



(22) 



(23) 



We note that I(x) plays the role of the width w(x) of a solid channel, that is, the integral 
in the previous equation becomes equal to the width of the channel if we model the rigid 
boundaries by a potential field that is zero (infinite) inside (outside) the channel. In fact, 
replacing I(x) by w(x) in Eq. fT22|) we obtain Eq. (25) in Ref. ll|. This suggests that 
the confining potential V(x, z) has similar effects on the transport of Brownian particles to 
those induced by a channel with solid walls. 

Before we proceed to solve the equation for fo(x) it is also interesting to compare it 
with the equation obtained from the direct extension of the Fick- Jacobs approximation in 
the presence of an external force. In fact, in the narrow geometry that we are considering 



here, we can assume that the particle will reach local equilibrium in the cross section fairly 
rapid compared to its diffusive or convective motion along the channel. This separation of 
time scales suggests that the conditional probability P(z/x, t) could be approximated by the 
equilibrium distribution P eq (z/x). Then, the total probability density takes the form, 



P(x, z, t) = P(z/x, t)P(x, t) « P eq (z/x)p(x, t) 
where p (x, t) is the marginal probability distribution, 

POO 

p(x,t) — / P(x,z,t)dz, 



-V(x,z) 



I(x) 



-p{x,t), 



(24) 



(25) 

J — oo 

Considering the long time limit, we can then use the same approximation for the asymptotic 
distribution, 

e -V(x,z) 

Poo(x,z) « -——p(x). (26) 
I(x) 

Finally, we can project the problem into the longitudinal direction by substituting this 
expression into the governing equation for the probability density and integrating over the 
cross-section, 



d_ 

dx 



Pe 



P 



Hx] 



d 



V 



dx \ I(x, 



Hx) 



0. 



(27) 



This is the form o: 



the Fick- Jacobs equation used in the presence of an external force 



in previous studies |5j, [9]. It is interesting to point out that, by making the substitution 
p(x) = fo(x)I(x), the previous equation becomes identical to Eq. (1221) . which was derived 
from the asymptotic analysis. Therefore, the leading order term of the asymptotic analysis 
is equivalent to the the proposed extension of the Fick- Jacobs approximation in the presence 
of an external field. Then, the validity of the ad-hoc position-dependent effective diffusivity 
proposed in Ref. [h]] could, in principle, be tested by extending the present approach to 
higher orders in the slenderness parameter. 

We now go back to equation (1221) to obtain the general solution for fo(x), 



^Pex 



Ci 



-Pe x 



Hx. 



-dx + Ci 



(28) 



where C\ and C2 are the constants of integration. The leading order contributions to the 
average velocity and the dispersion coefficient can then be calculated from the probability 
distribution. First, we obtain the general expression for the average velocity following a 
derivation analogous to that presented in Ref. 



10| (appendix A) for the particle current, 



(v) 



-Pe\ 



dxe Pex I(x) 



x+1 



dx 



3 Pe x 



I(x')\ 



(29) 



In order to calculate the dispersion coefficient given in Eq. (|T3|) we first propose a regular 
perturbation expansion for the -B-field, analogous to those proposed for the probability 
density and flux, 

B ~ B + e 2 B l + e 4 B 2 + ■■■ . (30) 

Substituting the expansion for both the probability density and the -B-field into Eq. ( jTTj) it 
is straightforward to show that the zeroth-order balance is an homogeneous equation that 
is satisfied by an arbitrary function B (x). The governing equation for B (x) is obtained 
from the 0(e 2 ) balance, after integrating over the cross-section and taking into account the 
no-flux boundary condition for the -B-field in the ^-direction, 

d*B 0+ f P ' -(v)\dB 1 ^ {v) (3l) 



dx 2 \ Po J dx 
The remaining boundary condition is, 

S (l) -So(0) = -1. (32) 
The above two equations determine Bq uniquely to within an arbitrary additive constant, 



which does not affect the calculation of the effective dispersion coefficient 



16| 



D j dx f po(x,z) dz. (33) 



TRANSPORT OF BROWNIAN PARTICLES CONFINED BY A PARABOLIC PO- 
TENTIAL 

In this section, we discuss a specific example of the confining potential that illustrates our 
previous results and allows direct comparison with numerical results. Consider the transport 
of Brownian particles confined to a narrow channel by a parabolic potential of periodically 
varying width, 

o \ 2 / \ 2 

2tt \ ( z — z 



V(x, z) = ^cos — x + dj y-J^J k sT, (34) 

where zq is the position of the center of the channel, which we assume to be constant (straight 
centerline), and d > 1 determines the minimum opening of the channel. The slenderness of 
the geometry is given by e 1, which is the ratio of the characteristic width of the potential 
in the ^-direction, eL, to the length of one period, L. Using the same dimensionless variables 
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introduced before and dropping again all bars for the sake of simplicity, we obtain, 

\2 



V(x,z) 



25 2 (x) 



where 



5(x) 



(35) 



(36) 



y/2 (cos2vra; + d) 

It is clear then that the particles will be confined in the z direction by a parabolic potential 
and that the width of the confining region is determined by 5(x). In fact, the equilibrium 
distribution of particles is given by the Boltzmann distribution with variance u(x) = 5(x). 
In Fig. Q] we plot the equipotential lines z — z = ±25 (x), with d = 1.2 and z = 8. In 
equilibrium, approximately 95% of the particles are confined to the region enclosed by these 
lines. 




FIG. 1: Equipotential lines corresponding to z = zq = ±25 (x) for a periodic potential V(x,z) = 
(cos 2ttx + d) 2 (z — zq) 2 with d = 1.2 and zq = 8. 

For this specific potential we obtain I(x) = \^2n5(x) = y^cos 2ttx ± d)~ x , and 

p / \ C\ ( 2ix sin 2irx — Pe cos 2irx d\ _ Ppr .„_. 
MX) = -V*{ ^UPe? ftj +C2e ' (37) 

Where C\ and C% are integration constants. In particular, C\ is the first constant of in- 
tegration of Eq. (|22l) . which corresponds to the total flux in the x-direction, C\ = (Jo). 
Thus, imposing periodicity and the normalization condition we obtain the leading order 
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probability density, 



(Jq) -v(x,z) ( 2vrsin27rx — Pecos27rx d 



Mx , z) = -^-W, _,, (38) 



and average velocity 



(Pe 2 + 4vr 2 ) Pe 
Pe 2 + 



W = (v) = \ 2 , J d • (39) 



Finally the leading order term of the marginal probability distribution, corresponding to 
the distribution in the heuristic extension to the Fick- Jacobs approximation, is given by 



p(x) 



Pe 2 



Pe 2 + 



4w 2 d 



1 27rPe sin 2wx — 47r 2 d 
Pe 2 d + cos(2vra) 



(40) 



Note that in the limiting case in which diffusive transport is dominant (vanishingly small 
Peclet number) we recover the equilibrium Boltzmann distribution, 



A2 _ 1 

lim p (x,z) = J e" v{x ' z \ (41) 

Pe^o V 7r 



On the other hand, in the limit of deterministic motion we obtain, 



v / \ d + cos(2irx) _ v(xz) ,. n . 
hm pq(x,z) = 7= -e v ( x ' z >. (42) 



In this case, there is no boundary layer developing due to the large magnitude of the driving 



force, as it would be the case in the presence of solid walls and permeating forces [15]. 
Therefore, the diffusive fluxes are negligible compared to convective transport, and the 
average velocity tends to its bulk value, something that is also clear from Eq. ( |39i) . In 



addition, the integral of the asymptotic distribution over a cross section is uniform, which is 
also clear from Eq. (|22|) . Note, however, tha 
satisfies e 2 Pe 1, as discussed in Ref. |ll| . 



also clear from Eq. (1221) . Note, however, that this limit is valid if the slenderness condition 



BROWNIAN DYNAMICS SIMULATIONS 

We performed Brownian dynamics simulations of particle transport under confinement 
by a potential landscape to examine the previous asymptotic results in more detail. The 
motion of Brownian particles in a viscous solvent in the limit of vanishingly small inertia is 
governed by the overdamped Langevin equations, 

dr f)V 

I ,_ = F __ + V ^ C(t) , ,43, 
10 



and 

where ((t) is a zero-mean Gaussian white noise with the correlation (Ot^iKjfe)) = 2SijS(ti — 
£2)- The value of the slenderness parameter used in the numerical simulations is e = 0.01. 




FIG. 2: Marginal probability density p(x) for different values of the Peclet number (Pe). The solid 
lines indicate the asymptotic results and the symbols correspond to the results of the Brownian 
dynamics simulations. 

In figure [2] we show the marginal probability distribution p(x) in a single period of the 
potential, for different values of the Peclet number Pe = 1, 10, 30, 50. The leading order term 
of the distribution, given by Eq. (flOl) . agrees well with the results of the Brownian dynamics 
simulations for all Peclet numbers. The figure also shows that the marginal probability 
distribution tends to a uniform distribution as the Peclet number increases. This indicates 
that, as discussed before, the effect of the potential on the particle distribution, as well as the 
contribution of diffusive transport to the average velocity of the particles, decreases as the 
Peclet number increases. Figure [3] shows that the average velocity of the particles obtained 
in the simulations is accurately described by the leading order term in the perturbation 
expansion. Finally, in figure H] we compare the effective dispersion coefficient obtained from 
the leading order term in the -B-field (numerically solving Eq. (|3~T|) by means of finite 
differences) with that computed directly from the simulations. We observe good agreement 
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FIG. 3: Average velocity along the channel as a function of the Peclet number. The solid line 
corresponds to the leading order term of the asymptotic results and the solid circles represent the 
values computed from the Brownian dynamics simulations. 

between the simulations and the asymptotic analysis. 




FIG. 4: Asymptotic dispersion coefficient as a function of the Peclet numbers. The solid line 
corresponds to the leading order term in the asymptotic analysis. The solid circles corresponds to 
the dispersion coefficient calculated from the Brownian dynamics simulations. 
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CONCLUSIONS 



We have investigated the transport of Brownian particles driven by a uniform external 
force. The particles are confined to a narrow periodic channel by a parabolic potential in 
the cross-section. We used asymptotic methods to obtain the leading order solution of the 
two-dimensional Smoluchowski equation for the long-time probability distribution of par- 
ticles reduced to a single period of the potential. We first showed that the leading order 
analysis reproduces a previously proposed extension of the Fick- Jacobs approximation to 
biased transport. We thus provide a systematic method to improve on the Fick- Jacobs ap- 
proximation through higher order analysis. We also showed that the leading order equation 
is equivalent to that obtained for solid channels and thus demonstrated that a confining po- 
tential has analogous effects on the distribution of particles and their transport parameters, 
such as the average velocity and the asymptotic dispersion coefficient. We then analyzed the 
case of a cosine variation in the aperture of the confining channel and compared the results 
with Brownian dynamics simulations. We compared the long-time marginal distribution as 
well as its first moments (average velocity and dispersion coefficient) and obtained excel- 
lent agreement with the first order in the asymptotic analysis over a wide range of Peclet 
numbers. 

This material is partially based upon work supported by the National Science Foundation 
under Grant No. CBET-0731032. 
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